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We revie-w the current status of strangeness as signature of the forma- 
tion and dissociation of the deconfined QGP at the SPS energy scale, and 
present the status of our considerations for RHIC energies. By analyzing, 
1-5 ! -within the frame-work of a Fermi statistical model, the hadron abundance 

and spectra, the properties of a disintegrating, hadron evaporating, decon- 

Q^l fined QGP fireball are determined and can be compared -with theory for 

D . the energy range 160-200AGeV on fixed target. We discuss in more detail 

i-Q ' our finding that the pion yields occur near to pion condensation condition. 

Dynamical models of chemical strangeness equilibration are developed and 
. . - applied to obtain strangeness production in a QGP phase at conditions 

KS ] found at SPS and expected at RHIC. The sudden QGP break up model 

that -works for the SPS data implies at RHIC dominance of both baryon, 
and antibaryon, abundances by the strange baryon and antibaryon yields. 

PACS numbers: 12.38.Mh, 25.75.-q, 25.75.Dw, 25.75.Ld 

1. Introduction 

Strange particle signatures for the formation and evolution of the de- 
confined quark-gluon phase of elementary matter (QGP) has been a subject 
developed quite intensely for the past 20 years. We revie-w here our progress 
since the last major revie-w [1], highlighting our analysis of the Pb-Pb data 
at SPS, and our predictions for hyperon yields from QGP at RHIC [2]. 

We first describe in section 2, the Fermi model [3] analysis of the multi- 
particle production processes in 158^ GeV Pb-Pb collisions carried out at 

^ Support by U.S. Department of Energy under grant DE-FG03-95ER40937 . 
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CERN-SPS. Strongly interacting particles are believed to be produced with 
a probability coniniensurating to the size of the accessible phase space. The 
numerical methods which have been developed in the context of an analysis 
of the lighter 200yl GeV S-Au/W/Pb system [4] are described. We have 
in particular shown [5, 6] that consideration of the light quark chemical 
non-equilibrium is necessary in order to arrive at a consistent interpretation 
of the experimental results of both the wide acceptance NA49-experiment 
[7, 8, 9, 10, 11] and central rapidity strange (multi) strange (anti)baryon 
WA97-experiment [12, 13, 14]. This resulted also in considerable reduction 
of the chemical freeze-out temperature: we find Tj = 145 it 5 MeV, while 
originally it has been estimated to be [15, 16] Tf = 180-290 MeV. 

Such a low freeze-out temperature is more consistent with the assump- 
tion we make that there is no change of hadronic particle abundance after 
the deconfined QGP source has dissociated. This scheme is called sudden 
hadronization [17, 18]. This can occur if hadronic particles are produced 
either in: 

a) an evaporation process from a hot expanding, surface or 

b) a sudden global hadronization process. 

Our sudden hadronization scheme works very well, and can be considered 
as established on view of many studies that could describe quite diverse 
data. Certain surprising features of Pb-Pb results that are seen within 
such analysis, and in particular the finding that the pion yield is governed 
by a fugacity that is close to the condensation point, as we shall show in 
section 3, lead us presently to favor the scenario b). That being the case 
one may further suppose that a super-cooled plasma occupies a relatively 
large spatial volume, and it undergoes a global explosive decomposition into 
individual hadrons, maximizing hadron occupancies and thus the entropy 
content in the confined phase at the near-pion-condensation condition. 

Pertinent results of our analysis of the Pb-Pb system are addressed in 
section 3, where we have reevaluated our current results in consideration of 
some small change of the experimental data. We address all available SPS 
NA49 and WA97 experimental data, except for and Q particles. It is im- 
portant to realize that if we succeed to describe well a particle yield within 
the Fermi model, it means that the majority of all particles of the partic- 
ular type is produced by the statistical mechanisms we address here. In 
principle there could be many other production mechanisms, and they add 
to the yields. Thus if our description fail, an acceptable failure is the one 
which under-predicts the yield. When statistical model predicts very little 
if any production, the other reaction pictures may indeed be dominant and 
we should at least hesitate in our attempt to describe all the rare particle 
yields. The prime candidate for such consideration and omission from sta- 
tistical analysis is the totally strange 0,{sss) and its antiparticle: they are 
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triply strangeness suppressed and are very heavy with Mq = 1672 MeV, thus 
again significantly suppressed, especially at low chemical freeze-out temper- 
ature. Their total statistical multiplicity is by a good distance smallest of 
all 'stable' hadrons. Consequently, their production pattern is easily altered 
by, e.g., in source strangeness clustering. We have found that if we use the 
results we obtain about the properties of the source [6] , in order to compute 
the yields of O and fi, we invariably see that we obtain only a fraction, 
40-50%, of all particles observed. Our preliminary conclusion is that we 
should NOT explore these particles in the statistical production model. 

Among results that we obtain in section 3 is the, on a first sight, surpris- 
ing overpopulation of the strangeness phase space occupancy. We explain 
how this can occur in section 4, where the kinetic theory for computation 
of the chemical strangeness flavor abundance equilibration is presented. We 
extend our past study of strangeness production at SPS conditions and show 
that, at the time of QGP breakup at RHIC energies, there is also in gen- 
eral full chemical equilibrium, indeed that one can expect over-saturation 
of strangeness flavor, just as at SPS. Our numerical study is based on the 
dynamics of the phase space occupancy rather than particle density, and we 
eliminate much of the dependence on the dynamical flow effects by incor- 
porating in the dynamics considered the hypothesis of entropy conserving 
matter flow and evolution. We will make two assumptions of relevance for 
the results we obtain: 

• the kinetic (momentum distribution) equilibrium is reached faster than 
the chemical (abundance) equilibrium [19, 20]; 

• gluons equilibrate chemically significantly faster than strangeness [21]. 
The first assumption allows to study only the chemical abundances, rather 
than the full momentum distribution, which simplifies greatly the structure 
of the master equations; the second assumption allows to focus after an 
initial time tq has passed on the evolution of strangeness population: tq is 
the time required for the development to near chemical equilibrium of the 
gluon population. As we shall see, the strange quark mass mg is the only 
undetermined parameter that enters strangeness yield calculations. The 
overpopulation of the strangeness phase space, seen in SPS data arises for 
relatively small ms(lGeV)~ 200 MeV. 

In the following section 5, we use the experience we have with the SPS 
systems and with the theoretical studies of strangeness production in QGP, 
in order to estimate the strange particle production that is likely to occur at 
RHIC. Some remarkable particle abundance results arise, since during the 
break-up of the QGP phase there is considerable advantage for strangeness 
fiavor to stick to baryons. This can be easily understood considering that 
production of strange baryons over kaons is favored by the energy balance, 
i.e. : E{K + vr) < i?(N+K). Since at RHIC most hadrons produced are 
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mesons, and baryons form just a small fraction of all particles, initially 
we expect and will show, in section 5, that hyperon production dominates 
baryon production, i.e., most baryons and antibaryons produced will be 
strange. A remarkable consequence of the sudden hadronization scenario 
is that this situation is maintained and thus hyperon dominance should be 
observed at RHIC. If indeed this prediction is born out in the experiment, 
it will prove that the there was formation of deconfined phase, followed by 
sudden hadronization. 

We note that at SPS energies described in section 3, there is still an 
appreciable relative baryon abundance among all hadrons (about 15%) and 
thus while hyperon dominance begins to set in, there are (literally speaking) 
still some non-strange baryons left. With increasing energy the yield of 
strange quark pairs per baryon rises, and at the same time the relative 
abundance of baryons among all hadrons diminishes, the relative population 
of non-strange baryons decreases rather rapidly and at RHIC energies the 
hyperons and/or antihyperons are the dominant strange particle fractions. 

We are not aware that other studies reported in literature about RHIC 
conditions have this remarkable result, see, e.g., [22]. It is thus interesting to 
record the two major quantitative differences of the behavior of deconfined 
matter we are considering: 

• in QGP the particle density is high enough to assure that the required 
abundance of strangeness can be actually produced [1, 23, 24, 25], while in 
hadron phase it was shown that, even at SPS energy, this is not the case [26]. 

• overpopulation of hadron phase space occupancies occurs naturally when 
the entropy rich QGP phase disintegrates into hadrons, which cannot be 
expected in hadron based kinetic reactions. 

2. Contemporary Fermi model of Hadron Production 

We use 6 parameters to characterize the spectra and abundances of par- 
ticles. Will describe these discussing their values, assuming a QGP source: 

1) The strange quark fugacity As = 1 can be obtained from the require- 
ment that strangeness balances [27]: 

{Ns-Ns) = 0. (1) 

However, the Coulomb distortion of the strange quark phase space plays an 
important role in the understanding of this constraint for Pb-Pb collisions 
[5], leading to the Coulomb-deformed value As = 1.10 , see also Eq. (10). 

2) Strange quark phase space occupancy 7s can be computed as we show 
in section 4 within the established kinetic theory framework for strangeness 
production [1, 23]. For a rapidly expanding system the production processes 
will lead to an oversaturated phase space with 7s > 1 . 
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3) The equilibrium phase space occupancy of hght quarks 7^ is expected 
to exceed unity significantly to accommodate the excess entropy content 
in the plasma phase [27]. There is an interesting constraint that arises if 
hadronization is sudden in the sense that particles are produced at the same 
time, forming for pions a Bose gas. As we shall discuss at the end of this 
section, see Eq. (9) this leads to an upper limit: 

7,<7^^^e™^/2^. (2) 

4) The collective surface expansion velocity should remain below the 
relativistic sound velocity [1]: 

Vc < I/VS. (3) 

5,6) If we assume that the stopping of the baryon number and energy 
is similar [1], we know the energy per baryon content in the reactions and 
then the equations of state produce a further constraint between chemical 
freeze-out temperature Tf and light quark fugacity Xq or equivalently, the 
baryochemical potential: 

^B = 3r/lnA,. (4) 

The difference between Aj and 7^ is that, e.g., for strange and anti- 
strange quarks the same factor 7^ applies, while the antiparticle fugacity is 
inverse of the particle fugacity. The proper statistical physics foundation 
of 7j is obtained considering the maximum entropy principle: it has been 
determined that while the limit 7j ^ 1 maximizes the specific chemical 
entropy, this maximum is extremely shallow, indicating that a system with 
dynamically evolving volume will in general find more effective paths to 
increase entropy, than offered by the establishment of the absolute chemical 
equilibrium [28]. 

The abundances of the final state particles is most conveniently described 
by considering the phase space distribution of particles. The relative number 
of primary particles freezing out from a source is obtained noting that the 
fugacity and phase space occupancy of a composite hadronic particle is ex- 
pressed by its constituents and that the probability to find all j-components 
contained within the i-th emitted particle is: 

Taking the Laplace transform, we find, e.g., for the strange sector, to 
the following partition function, like expression: 

InZ, = ^-^ {(A,A-i + X;'Xghs^gCl,FK + (A.A^ + Xj' X-^hs^^ClFy 
+{X',Xq + Xj^Xq'hh,CI,Fs + {Xl + Xj^h^,CI,Fn} , (6) 
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where the kaon, hyperon, cascade and omega degrees of freedom are in- 
cluded. The phase space factors Fj of the strange particles are (with gi 
describing the statistical degeneracy): 



F^ = Y.^'^^'^^J^)- (7) 



In the resonance sums X^i ^ known strange hadrons should be counted. 
The function W{x) = x'^K2{x), where K2 is the modified Bessel function, 
arises from the phase-space integral of the different particle distributions 
fip). It is important to remember that this expression does not describe 
the properties of a gas of hadrons, thus it is not a partition function, even 
if we give the Laplace transform of the phase space such formal semblance. 
When the source of the particles is subject to flow, the Laplace transform 
that leads to the above expression is considerably more involved, the spectra 
and thus also multiplicities of particles emitted are obtained replacing the 
Boltzmann factor in Eq. (5) by [29]: 



e 



-^'/^ - ^Jdn,^,{l-v,-p,/E,)e''^('-'^-P^^/^'\ 



1 



7c = ^r^f, (8) 



-^c 



a result which can be intuitively obtained by a Lorentz transformation be- 
tween an observer on the surface of the fireball, and one at rest in laboratory 
frame. In certain details the results we obtain confirm the applicability of 
this simple approach. We consider for SPS energy range the radial flow 
model, perhaps the simplest of the complex flow cases possible , but it 
suffices to fully assess the impact of flow on our analysis. 

While the integral over the entire phase space of the flow spectrum yields 
as many particles with and without flow, when acceptance cuts are present 
particles of different mass experience differing flow effects. Here, we note 
that the final particle abundances measured in an experiment are obtained 
after all unstable hadronic resonances are allowed to disintegrate and feed 
the stable hadron spectra. In order to minimize the impact of unknown flow 
pattern at hadron freeze-out when considering particle abundances mea- 
sured in a restricted phase space domain, we study particle abundance ra- 
tios involving what we call compatible hadrons: these are particles likely to 
be impacted in a similar fashion by collective flow dynamics in the fireball. 

We now return to review the case of pions, which is exceptional since we 
will be considering a rather large values of 7^ > 1.5. The chemical fugacity 
for a particle composed of a light quark-antiquark pair is 7^ . Thus the Bose 
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142 [MeV] 



S/V [fm-3] 
N/V [fm-3] 
E/V [GeV fm-3] 




Fig. 1. Dependence of pion gas properties A^/F-particle, E/V-eneigy and S/V- 
entropy density, as function of 7^ at T = 142 MeV. 



distribution in momentum space has the shape: 



UE) 



^-2eE./r 



1 



E^ 



m^ +p2 



(9) 



The range of values for 7^ is bounded from above by the Bose singularity. 
When 7g -^ 7^, see Eq. (2), the lowest energy state (in the continuum limit 
with p ^ ) will acquire macroscopic occupation and a pion condensate is 
formed. Such a condensate 'consumes' energy without consuming entropy of 
the primordial high entropy QGP phase. Thus a condensate is not likely to 
develop, but the hadronization process may have the tendency to approach 
the limiting value in order to more efficiently connect the deconfined and 
the confined phases, since, as we show in Fig. 1, the entropy density is nearly 
twice as high at 7^ ~ 7^ than at 7^ = 1. 

To see clearly how this can occur, we looked more closely at the relative 
properties of a pion gas for 7<j ^ 7^ • In Fig. 2, we see the relative change in 
energy per pion, (inverse of) entropy per pion, and energy per entropy, for 
fixed T = 142 MeV corresponding to our best fit condition. We see that a 
hadronizing gas will consume at higher 7^ less energy per particle, and that 
the energy per entropy is nearly constant. Dissociation into pions at 7^ -^ 7^ 
appears thus to be an effective way to convert excess of entropy in the plasma 
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Fig. 2. Dependence of pion gas properties (iV-particle, S-energy and S'-entropy) 
density as function of 7^ for T = 142 MeV. a) ratios relative to equilibrium value 
7q = 1; b) relative ratios, thus E/N, S/N and E/S. 



into hadrons, without need for reheating, or a mixed phase which would 
allow the volume to grow. In short, the finding of the maximum allowable 
7q is intrinsically consistent with the notion of an explosively disintegrating 
QGP phase. 



3. Update of SPS Experimental Data Analysis 

In the past year our work addressed our discovery that consideration 
of the light quark chemical non-equilibrium is necessary in order to arrive 
at a consistent interpretation of the experimental results emanating from 
CERN [30]. We have also incorporated in our earlier analysis of the Pb-Pb 
system [5] a study of collective matter flow. Properties of the dense fireball 
as determined in this approach offer clear evidence that a QGP disintegrates 
at Tf ~ 144 MeV, corresponding to energy density e = 0(0.5) GeV/fm'^ [31]. 
With flow, the analysis addresses also the m^-slopes of strange particles. 
Notably, the near equality of (inverse) slopes of nearly all strange baryons 
and antibaryons arises by means of the sudden hadronization at the surface 
of an exploding QGP fireball. In the hadron based microscopic simulations 
this behavior of m^-slopes can also arise allowing for particle-dependent 
freeze-out times [33]. 

We note that though we use six parameters to characterize the hadron 
phase space at chemical freeze-out, compare section 2, there are only two 
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truly unknown properties: the chemical freeze-out temperature Tf and the 
light quark fugacity Ag (or equivalently, the baryochemical potential Eq. (4)) 
— we recall that the parameters 7j, i = q,s controls overall abundance of 
quark pairs, while Aj controls the difference between quarks and anti-quarks 
of given flavor. As already noted earlier, the four other parameters are not 
arbitrary, and we could have used their tacit and/or computed values: 

1) the strange quark fugacity A^ is usually fixed by the requirement that 
strangeness balances {s — s) = [18]. The Coulomb distortion of the strange 
quark phase space plays an important role in the understanding of this con- 
straint for Pb-Pb collisions, see Eq. (10) [5]; 

2) the strange quark phase space occupancy 7^ can be computed within 
the established kinetic theory framework for strangeness production [1, 23]; 

3) the tacitly assumed equilibrium phase space occupancy of light quarks 

7g = 1; 

4) assumed collective expansion to proceed at the relativistic sound ve- 
locity, Vc = l/\/3 [1]. 

However, the rich particle data basis allows us to find from experiment the 
actual values of these four parameters, allowing to confront the theoretical 
results and/or hypothesis with experiment. 

The value of A^ we obtain from the strangeness conservation condition 
(s — s) = in QGP is, to a very good approximation [5]: 



A,^aX^3 = 1, Aq^ -'^- ^, . (10) 



^ " ' "^^ k,d' 



r 



Aq < 1 expresses the Coulomb deformation of strange quark phase space. 
This effect is relevant in central Pb-Pb interactions, but not in S-Au/W/Pb 
reactions. Aq is not a fugacity that can be adjusted to satisfy a chemical 
condition, since consideration of Aj, i = u,d,s exhausts all available chem- 
ical balance conditions for the abundances of hadronic particles. The sub- 
script Rf in Eq. (10) reminds us that the classically allowed region within 
the dense matter fireball is included in the integration over the level density. 
Choosing R{ = 8fm,T = 140 MeV, m^ = 200 MeV (value of 7^ is practically 
irrelevant), for Zf = 150 the value is A^ = 1.10 . 

The available compatible particle yield ratios (excluding 0, and 0,, see 
section 1) are listed in table 1, top section from the experiment WA97 for 
p± > 0.7 GeV within a narrow Ay = 0.5 central rapidity window. Fur- 
ther below are shown results from the large acceptance experiment NA49, 
extrapolated to full 47r phase space coverage. We first fit 11 experimental 
results shown in table 1, and than turn to include also the m_L-slope in our 
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Table 1. WA97 (top) and NA49 (bottom) Pb-Pb 158^ GeV particle ratios and 
some of our theoretical results, see text for explanation. 



Ratios 


Rcf. 


Exp. Data 


Pbjo 


Pbj^ 


pi^lsb 


pi^lsc 


H/A 


13 


0.099 ± 0.008 


0.104 


0.103 


0.105 


0.103 


S/A 


fl3l 


0.203 ± 0.024 


0.214 


0.208 


0.209 


0.206 


A/A 


[131 


0.124 ± 0.013 


0.124 


0.125 


0.124 


0.125 


" /" 


[13] 


0.255 ± 0.025 


0.256 


0.252 


0.248 


0.251 


(5+5) 
(A+A) 


[34] 


0.13 ± 0.03 


0.126 


0.122 


0.124 


0.122 


KV<t^ 


|8| 


11.9 ± 1.5 


14.2 


13.3 


13.0 


13.4 


K+/K^ 


[9] 


1.80± 0.10 


1.80 


1.82 


1.78 


1.83 


P/P 


[7] 


18.1 ±4. 


17.3 


16.7 


16.6 


16.6 


A/p 


[35] 


3. ± 1. 


2.68 


2.11 


2.11 


2.11 


K^/B 


[361 


0.183 ± 0.027 


0.181 


0.181 


0.163 


0.188 


h-/B 


[371 


1.97 ± 0.1 


1.96 


1.97 


1.97 


1.96 




X'r 




3.6 


2.5 


3.2 


2.6 




N;p;r 




11;5;2 


12;6;2 


12;5;2 


12;5;2 



considerations, and thus have 12 data points. The total error: 



Xt 



E,« 



^ ''cxp J 



(Ai?; 



:ii) 



cxp; 



for the four theoretical columns is shown at the bottom of this table along 
with the number of data points W, parameters 'p' used and (algebraic) 
redundancies 'r' connecting the experimental results. For r 7^ it is more 
appropriate to quote the total Xt' with a initial qualitative statistical rel- 
evance condition Xt/{^ ~ p) < 1- The first theoretical columns refer to 
results without collective velocity Vc (subscript 0) the three other with fit- 
ted Vc (subscript Vc)- In column three, superscript 's6' means that A^ is 
fixed by strangeness balance and, in column four, superscript ^sc^ means 
that 7g = 7q = e'^^''^'^f, that is 7^ is fixed by its upper limit, the pion 
condensation point. All results have been newly recomputed, to account for 
slightly higher value of the ratio h^/B [37]. 

It is interesting to note that the highest confidence result is obtained in 
the last column, just when the light quark phase space occupancy assumes 
value at the pion condensation point: here the number of degrees of free- 
dom is higher than in the second column, obtained without constraint. It is 
unclear at present what is the full extent of this remarkable result. Another 
interesting insight is that radial flow always on its own improves our abil- 
ity to describe the data. However, m± spectra offer another independent 
measure of flow, and confirm very strongly our findings about the value of 
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Table 2. Experimental and theoretical m± spectra inverse slopes Tth- Left Pb-Pb 
results from experiment NA49 [11] for kaons and from experiment WA97 [14] for 
baryons; right S-W results from WA85 [38]. 





T^^ MeV; 


Tt? MeV 


Tl MeV 


^t'h MeV 


rpK-> 


223 ± 13 


241 


219 ± 5 


215 


rpA 


291 ± 18 


280 


233 ± 3 


236 


rpA 


280 ± 20 


280 


232 ± 7 


236 


T= 


289 ± 12 


298 


244 ± 12 


246 


T" 


269 ± 22 


298 


238 ± 16 


246 



Vc- We proceeded as follows: for a given pair of values Tj and Vc we eval- 
uate the resulting mj_ particle spectrum and analyze it using the spectral 
shape and kinematic cuts employed by the experimental groups. To find the 
best values we consider just one 'mean' strange baryon experimental value 
T^ = 260 lb 10, since within the error the high m± strange (anti)baryon 
inverse slopes are overlapping. Thus when considering Vc along with Tj_ we 
have one parameter and one data point more. Once we find best values of 
Tf and Vc, we study again the inverse slopes of individual particle spectra. 
We obtain an acceptable agreement with the experimental T| as shown in 
left section of table 2 . 

For comparison, we have also considered in the same framework the 
S-induced reactions, and the right section of table 2 shows a good agree- 
ment with the WA85 experimental data [38]. We used as the 'mean' ex- 
perimental slope data point T^ = 235 it 10. We can see that within a 
significantly smaller error bar, we obtained an accurate description of the 
rrij^-slope data. This analysis implies that the kinetic freeze-out, where elas- 
tic particle-particle collisions cease, cannot be occurring at a condition very 
different from the chemical freeze-out. However, one pion HBT analysis at 
p± < 0.5 GeV suggests kinetic pion freeze-out at about Tk ~ 120 MeV [39]. 
A possible explanation of why here considered p± > 0.7 GeV particles are 
not subject to a greater spectral deformation after chemical freeze-out, is 
that they escape before the bulk of softer hadronic particles is formed. At 
least for strange brayons and antibaryons this is the result also seen in a 
recent microscopic study of the freeze-out process [22] . 

The six statistical parameters describing the particle abundances are 
shown in the top section of table 3, where we also show in the last column 
for comparison, the best result for S-induced reactions, where the target has 
been W/Au/Pb [4]. The errors in the results are one standard deviation 
errors arising from the propagation of the experimental measurement error, 
but apply only when the theoretical model describes the data well. All 
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Table 3. Heading: Xt^ number of data points N, parameters p and redundancies r; 
upper section: statistical model parameters which best describe the experimental 
results for Pb~Pb data, and in last column for S-Au/W/Pb data presented in 
Ref. [4] . Bottom section: specific energy, entropy, anti-strangeness, net strangeness 
of the full hadron phase space characterized by these statistical parameters. In 
column three we fix As by requirement of strangeness conservation, and in column 
four we choose jq = 7^, the pion condensation point. 







Pb|. 


Pb|f 


Pb|- 


S|. 


x^; N 


;p;r 


2.5; 12; 6; 2 


3.2; 12; 5; 2 


2.6; 12; 5; 2 


6.2; 16; 6; 6 


Tf MeV 


142 ± 3 


144 ± 2 


142 ± 2 


144 ± 2 


Vc 




0.54 ± 0.04 


0.54 ± 0.025 


0.54 ± 0.025 


0.49 ± 0.02 


\ 




1.61 ± 0.02 


1.605 ± 0.025 


1.615 ± 0.025 


1.51 ± 0.02 


As 




1.09 ± 0.02 


1.10* 


1.09 ± 0.02 


1.00 ± 0.02 


Iq 




1.7 ± 0.5 


1.8± 0.2 


Y = f,ra^/2Tf 


1.41 ± 0.08 


Ishq 




0.79 ± 0.05 


0.80 ± 0.05 


0.79 ± 0.05 


0.69 ± 0.03 


Ef/B 




7.8 ± 0.5 


7.7 ± 0.5 


7.8 ± 0.5 


8.2 ± 0.5 


Sf/B 




42 ± 3 


41 ± 3 


43 ± 3 


44 ± 3 


Sf/B 




0.69 ± 0.04 


0.67 ± 0.05 


0.70 ± 0.05 


0.73 ± 0.05 


{Sf-S 


f)/B 


0.03 ± 0.04 


0* 


0.04 ± 0.05 


0.17 ± 0.05 



results shown in table 3 have convincing statistical confidence level. For 
the S-induced reactions the number of redundancies r shown in heading of 
the table 3 is large, since same data comprising different kinematic cuts has 
been included in the analysis. It is quite reassuring that within error the 
freeze-out temperature Tf seen in table 3, is the same for both the S- and 
Pb-induced reactions, even though the chemical phase space occupancies 
differ greatly. This must be the case within our model of sudden freeze- 
out and constitutes its firm confirmation. The variation in the shape of 
the particle spectra is fully explained by a change in the collective velocity, 
which rises from v^ = 0.49 ± 0.02 to v^^ = 0.54 ± 0.04 ~ l/^/2, = 0.577. 
The light quark fugacity \q implies that baryochemical potential is /U^^ = 
203 lb 5 > //^ = 178 lb 5 MeV. As in S-induced reactions where A^ = 1, 
now in Pb-induced reactions, a value A^ — 1.1 characteristic for a source 
of freely movable strange quarks with balancing strangeness, i.e.. As = 1 is 
obtained, see Eq. (10). 

The values of 7^ > 1, seen in table 3, imply that there is phase space 
over-abundance of light quarks, to which, e.^., gluon fragmentation at QGP 
breakup prior to hadron formation contributes. 7^ assumes in our data anal- 
ysis a value near to where pions could begin to condense [30], Eq. (2). We 
found studying the ratio h~ /B separately from other experimental results 
that the value of 7^ ~ 7^ is fixed consistently and independently both, by 
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the negative hadron {h^), and the strange hadron yields. The unphysical 
range 7g > 7^ ~ 1.63 can arise (see column Pb|^^) since, up to this point, 
we use only a first quantum (Bose/Fermi) correction. However, when Bose 
distribution for pious is implemented, which requires the constraint 7^ < 7^, 
we obtain practically the same results, as shown in second column of table 3. 
We then show in table 3 the ratio "fs/lq — 0.8, which corresponds (approxi- 
mately) to the parameter 7^ when 7^ = 1 had been assumed. We note that 
7^^ > 1. This strangeness over-saturation effect could arise from the effect 
of gluon fragmentation combined with early chemical equilibration in QGP, 
7s (t < tf) ~ 1. The ensuing rapid expansion preserves this high strangeness 
yield, and thus we find the result 7^ > 1 , as is shown in Fig. 33 in [1]. 

We show, in the bottom section of table 3, the energy and entropy 
content per baryon, and specific anti-strangeness content, along with specific 
strangeness asymmetry of the hadronic particles emitted. The energy per 
baryon seen in the emitted hadrons is nearly equal to the available specific 
energy of the colhsion (8.6 GeV for Pb-Pb, 8.8-9 GeV for S-Au/W/Pb). 
This implies that the fraction of energy deposited in the central fireball must 
be nearly the same as the fraction of baryon number. The small reduction of 
the specific entropy in Pb~Pb compared to the lighter S-Au/W/Pb system 
maybe driven by the greater baryon stopping in the larger system, also seen 
in the smaller energy per baryon content. Both collision systems freeze out 
at energy per unit of entropy E/S = 0.185 GeV. There is a loose relation 
of this universality in the chemical freeze-out condition with the suggestion 
made recently that particle freeze-out occurs at a fixed energy per baryon 
for all physical systems [40], since the entropy content is related to particle 
multiplicity. The overall high specific entropy content we find agrees well 
with the entropy content evaluation made earlier [27] for the S-W case. 

Inspecting Fig. 38 in [1], we see that the specific yield of strangeness we 
expect from the kinetic theory in QGP is at the level of 0.75 per baryon, in 
agreement with the results of present analysis shown in table 3. This high 
strangeness yield leads to the enhancement of multi-strange (anti)baryons, 
which are viewed as important hadronic signals of QGP phenomena [41], 
and a series of recent experimental analysis has carefully demonstrated com- 
paring p-A with A-A results that there is quite significant enhancement 
[14, 42], as has also been noted before by the experiment NA35 [43]. The 
strangeness imbalance seen in the asymmetrical S-Au/W/Pb system (bot- 
tom of table 3) could be a real effect arising from hadron phase space prop- 
erties. However, this result also reminds us that though the statistical errors 
are very small, there could be in this asymmetric system a considerable sys- 
tematic error due to presence of a significant spectator matter component. 
In the symmetric Pb-Pb collisions this effect disappears, despite the fact 
that the freeze-out flow pattern could be much more complex and there 
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could be a distortion of particle spectra at low momenta not accounted for 
in our study, for we did not model the ratio of kaons to hyperons. Consid- 
ering this limitation it is indeed remarkable, how well the conservation of 
strangeness condition is satisfied, when it is not being enforced. 

4. RHIC and Dynamics of Strangeness Production 

In some key aspects, the methods to describe strangeness production 
which we have been developing diff'er from those obtained in other studies 
of chemical equilibration of quark flavor, in particular for RHIC conditions 
[44, 45, 46]. For example, we use running QCD parameters (both coupling 
and strange quark mass) with strong coupling constant Ug as determined at 
the M^o energy scale. We also incorporate entropy conserving flow into the 
dynamical equations directly, exploiting significant cancellations that occur, 
and thus obtain a relatively simple dynamical model for the evolution of the 
phase space occupancy 7^ of strange quarks in the expanding QGP. 

The phase space distribution fg can be characterized by a local temper- 
ature r(x, t) of a (Boltzmann) equilibrium distribution f^ , with normal- 
ization set by a phase space occupancy factor: 

fMx;t))^^s{T)fTip-T). (12) 

Eq. (12) invokes in the momentum independence of 7^ our first assumption. 
More generally, the factor 7j, i = g,q,s,c, allows a local density of gluons, 
light quarks, strange quarks and charmed quarks, respectively not to be de- 
termined by the local momentum shape, but to evolve independently. With 
variables (t, x) referring to an observer in the laboratory frame, the chemical 
evolution can be described by the strange quark current non-conservation 
arising from strange quark pair production described by a Boltzmann colli- 
sion term: 

+ pS)p,{t){avYT^^''-ps{t)p,{t) (a^)^^^^^'«^.(13) 

The factor 1/2 avoids double counting of gluon pairs. The implicit sums 
over spin, color and any other discreet quantum numbers are combined 
in the particle density p = J2s c ...I (^'^P f^ ^^^ ^^ have also introduced 
the momentum averaged production/annihilation thermal reactivities (also 
called 'rate coefficients'): 

, \ _ I d^Pi I d^P2cri2Vi2f{pi,T)f{p2,T) 

^'"''"''^^ I d^Pi I d^P2m, T)fiP2, T) ■ ^'^> 
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f{pi,T) are the relativistic Boltzmann/Jiittner distributions of two colliding 
particles i = 1,2 of momentum pi. 

The current conservation used above in the laboratory 'Eulerian' formu- 
lation can also be written with reference to the individual particle dynamics 
in the so called 'Lagrangian' description: consider ps as the inverse of the 
small volume available to each particle. Such a volume is defined in the local 
frame of reference for which the local flow vector vanishes 'y(^)^)llocal ~ ^' 
The considered volume 6Vi being occupied by small number of particles 6N 
{e.g., 5N = 1), we have: 

5Ns = PsSVi . (15) 

The left hand side (LHS) of Eq. (13) can be now written as: 

dps dvps _ 1 d6Ns _ dps I d5Vi 

dt dx ^ 6Vi dt ~ dt ^'5Vi dt ' ^ ' 

Since 6N and 6Vidt are L(orentz)-invariant, the actual choice of the frame 
of reference in which the right hand side (RHS) of Eq. (16) is studied is 
irrelevant and we drop henceforth the subscript /. 

We can further adapt Eq. (16) to the dynamics we pursue: we introduce 
p^{T) as the (local) chemical equilibrium abundance of strange quarks, 
thus p = jsPf- We evaluate the equilibrium abundance 6N^ = 5Vp'^{T) 
integrating the Boltzmann distribution: 

SNr' = [6VT^]^z'K,{z), z = -^, (17) 

where Ki, is the modified Bessel function of order u; we will below use: 
d[z''K^{z)]/dz = -z^K^-i . The first factor on the RHS in Eq. (17) is a 
constant in time should the evolution of matter after the initial pre-thermal 
time period tq be entropy conserving [47], and thus 5VT^ = (5Vo^o =Const. . 
We now substitute in Eq. (16) and obtain 

where T = dT/dt. Note that in Eq. (18) only a part of the usual flow- 
dilution term is left, since we implemented the adiabatic volume expansion, 
and study the evolution of the phase space occupancy in lieu of particle 
density. The dynamics of the local temperature is the only quantity we 
need to model. 

We now return to study the collision terms seen on the RHS of Eq. (13). 
A related quantity is the (L-invariant) production rate ^^12^34 q£ particles 
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per unit time and space, defined usually with respect to chemically equili- 
brated distributions: 

1 



A 



12^34 



1 + 6 



1,2 



-pTpf{^sVl2 



,12^34 



(19) 



The factor 1/(1 + ^1^2) is introduced to compensate double-counting of iden- 
tical particle pairs. In terms of the L-invariant A , Eq. (13) takes the form: 



Tp° 






+l<i{r)lq{T)A 



qq—fss_ 



jsir)Tsir)iA''-''+A' 



SS— >(jr(j'^ 



.(20) 



Only weak interactions convert quark flavors, thus, on hadronic time scale, 
we have 7s,g(T) = 75,g(T). Moreover, detailed balance, arising from the time 
reversal symmetry of the microscopic reactions, assures that the invariant 
rates for forward/backward reactions are the same, specifically 

^12^34 _ ^34^12 
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and thus: 



Tp', 



(dru 

\dT ' T~K2{z) 






72 (r)^^^--" 



7?(r) 



+^l{T)A'i^^'' 



1 



7.2 (r) 



(22) 



When all 7j ^ 1, the Boltzmann collision term vanishes, we have reached 
equilibrium. 

As discussed, the gluon chemical equilibrium is thought to be reached at 
high temperatures well before the strangeness equilibrates chemically, and 
thus we assume this in what follows, and the initial conditions we will study 
refer to the time at which gluons are chemically equilibrated. Setting A^ = 1 
(and without a significant further consequence for what follows, since gluons 
dominate the production rate, also Xq = 1) we obtain after a straightforward 
manipulation the dynamical equation describing the evolution of the local 
phase space occupancy of strangeness: 



2t,T 



t^ 7s Kijz) 
dT T^K2{z) 



1 



7s 



(23) 



Here, we defined the relaxation time Tg of chemical (strangeness) equilibra- 
tion as the ratio of the equilibrium density that is being approached, with 
the rate at which this occurs: 

^ -" (24) 



2 (A 



gg^ss 



Pf 

+ Ai^- 



+ 



ss993 printed on February 1, 2008 17 



The factor 1/2 is introduced by convention in order for the quantity r^ to 
describe the exponential approach to equihbrium. 

Eq. (23) is our final analytical result describing the evolution of phase 
space occupancy. Since one generally expects that 7^ ^ 1 in a monotonic 
fashion as function of time, it is important to appreciate that this equa- 
tion allows 7s > 1: when T drops below ?n,s, and l/r^ becomes small, the 
dilution term (2nd term on LHS) in Eq. (23) dominates the evolution of 
7s . In simple terms, the high abundance of strangeness produced at high 
temperature over-populates the available phase space at lower temperature, 
when the equilibration rate cannot keep up with the expansion cooling. This 
behavior of 7^ has been shown in [32, Fig. 2] for the SPS conditions with 
fast transverse expansion. Since we assume that the dynamics of transverse 
expansion of QGP is similar at RHIC as at SPS, we will obtain a rather 
similar behavior for 7^ . We note that yet a faster transverse expansion than 
considered here could enhance the chemical strangeness anomaly. 

Ts{T) , Eq. (24), has been evaluated using pQCD cross section and em- 
ploying NLO (next to leading order) running of both the strange quark 
mass and QCD-coupling constant q^ [48]. We believe that this method pro- 
duces a result for q^ that can be trusted down to 1 GeV energy scale which 
is here relevant. We employ results obtained with ^^(M^o) = 0.118 and 
ms(lGeV) = 200 MeV; we have shown results with ms(lGeV) = 220 MeV 
earlier [2]. There is some systematic uncertainty due to the appearance of 
the strange quark mass as a fixed rather than running value in both, the 
chemical equilibrium density pf^ in Eq. (24), and in the dilution term in 
Eq. (23). We use the value ms(l GeV), with the 1 GeV energy scale chosen 
to correspond to typical interaction scale in the QGP at temperatures under 
consideration. 



5. Expectations for Strange Hadron Production at RHIC 

We now combine our advances in theoretical models of strangeness pro- 
duction and data interpretation at SPS energies with the objective of mak- 
ing reliable predictions for the RHIC energy range [2]. First we address the 
question how much strangeness can be expected at RHIC. Numerical study 
of Eq. (23) becomes possible as soon as we define the temporal evolution of 
the temperature for RHIC conditions. We expect that a global cylindrical 
expansion should describe the dynamics: aside of the longitudinal flow, we 
allow the cylinder surface to expand given the internal thermal pressure. 
SPS experience suggests that the transverse matter flow will not exceed the 
sound velocity of relativistic matter v^ ~ c/v3. We recall that for pure 
longitudinal expansion local entropy density scales as S" oc T^ oc l/r, [47]. 
It is likely that the transverse flow of matter will accelerate the drop in en- 
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tropy density. We thus consider the fohowing temporal evolution function 
of the temperature: 



T{t) = To 



1 



[l+T2c/d){l+TV^/Rl_Y 



1/3 



(25) 



We take the thickness of the initial collision region at Tq = 0.5 GeV to be 
(i(To = 0.5)/2 = 0.75 fm, and the transverse dimension in nearly central 
Au-Au collisions to be i?_L = 4.5 fm. The time at which thermal initial 
conditions are reached is assumed to be tq = Ifm/c. When we vary Tq, the 
temperature at which the gluon equilibrium is reached, we also scale the 
longitudinal dimension according to: 

d(To) = (0.5 GeV/To)3l.5 fm . (26) 

This assures that when comparing the different evolutions of 7s we are 
looking at an initial system that has the same entropy content by adjusting 
its initial volume Vq. The reason we vary the initial temperature Tq down to 
300 MeV, maintaining the initial entropy content is to understand how the 
assumption about the chemical equilibrium of gluons, reached by definition 
at To, impacts our result. In fact when considering decreasing Tq (and thus 
increasing Vq), what we are doing is to begin the thermal production at a 
later time in the history of the collision. 

The numerical integration of Eq. (23) is started at tq, and a range of ini- 
tial temperatures 300 < Tq < 600, varying in steps of 50 MeV. The high limit 
of the temperature we explore exceeds somewhat the 'hot glue scenario' [19], 
while the lower limit of Tq corresponds to the more conservative estimates of 
possible initial conditions [47] . Since the initial p-p collisions also produce 
strangeness, we take as an estimate of initial abundance a common initial 
value 7s(To) = 0.15. The time evolution in the plasma phase is followed up 
to the break-up of QGP. This condition we establish in view of our analysis 
of the SPS results. We recall that SPS-analysis showed that the system 
dependent baryon and antibaryon m_|_-slopes of particle spectra are result 
of differences in collective flow in the deconfined QGP source at freeze-out. 
There is a universality of physical properties of hadron chemical freeze-out 
between different SPS systems, and in our analysis a practical coincidence 
of the kinetic freeze-out conditions with the chemical freeze-out. We thus 
expect, extrapolating the phase boundary curve to the small baryochemical 
potentials, that the QGP break-up temperature T? ~ 145 it 5 MeV will 
see just a minor upward change to the value T? ~ 150 it 5 MeV. 

With the freeze-out condition fixed, one would think that the major 
uncertainty in our approach comes from the initial gluon equilibration tem- 
perature Tq, and we now study how different values of Tq influence the final 
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().;! 0.4 

T [GeV] 

Fig. 3. Evolution of QGP-phase strangeness phase space occupancy 7s. a) as 
function of time and, b) as function of temperature for ?Tis(l GeV) = 200 MeV, see 
text for details. 



state phase space occupancy. We integrate numerically Eq. (23) and present 
7s as function of both time t in Fig. 3a, and temperature T in Fig. 3b, up 
to the expected QGP breakup at T^Hic ~ 150 ± 5 MeV. We see that: 

• widely different initial conditions (with similar initial entropy content) lead 
to rather similar chemical conditions at chemical freeze-out of strangeness, 

• despite a series of conservative assumptions, we find, not only, that strange- 
ness equilibrates, but indeed that the dilution effect allows an overpopula- 
tion of the strange quark phase space. For a wide range of initial condi- 
tions, we obtain a narrow band 1.15 > 7s(Tj) > 1 . We will in the follow- 
ing, taking into account some contribution from hadronization of gluons in 
strange/antistrange quarks, adopt what the value ^s{Tf) = 1.25. 

We now consider how this relatively large value of 7^, characteristic 
for the underlying QGP formation and evolution of strangeness, impacts 
the strange baryon and anti-baryon observable emerging in hadronization. 
Remembering that major changes compared to SPS should occur in rapid- 
ity spectra of mesons, baryons and antibaryons, we will apply the same 
hadronization model that worked in the analysis of the SPS data. This 
hypothesis can be falsified easily, since we expect, based and compared to 
the Pb-Pb 158^ GeV results: _ _ 

a) shape identity of all RHIC m± and y spectra of antibaryons p , A , S , 
since in our approach there is no difference in their production mechanism, 
and the form of the spectra is determined in a similar way by the local 
temperature and flow velocity vector; 
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b) the mj_-slopes of these antibaryons should be very similar to the 
result we have from Pb-Pb 158 A GeV since only a slight increase in the 
freeze-out temperature occurs, and no increase in collective transverse flow 
is expected. 

The abundances of particles produced from QGP within the sudden 
freeze-out model are controlled by several parameters we addressed ear- 
lier: the light quark fugacity 1 < Ag < 1.1 , value is limited by the expected 
small ratio between baryons and mesons (baryon-poor plasma) when the en- 
ergy per baryon is above 100 GeV, strangeness fugacity As ~ 1 which value 
for locally neutral plasma assures that {s — s) = 0; the light quark phase 
space occupancy 7g ~ 1.5, overabundance value due to gluon fragmentation. 
Given these narrow ranges of chemical parameters and the freeze-out tem- 
perature Tf = 150 MeV, we compute the expected particle production at 
break-up. In general we cannot expect that the absolute numbers of parti- 
cles we find are correct, as we have not modeled the important effect of flow 
in the laboratory frame of reference. However, ratios of hadrons subject to 
similar flow effects (compatible hadrons) can be independent of the detailed 
final state dynamics, as the results seen at SPS suggest, and we will look at 
such ratios more closely. 

Taking 7^ = 1.5 ^0.25 we choose the value of Ag, see the header of table 4, 
for which the energy per baryon (E/B) is similar to the collision condition 
(100 GeV), which leads to the range Ag = 1.03ib0.005. We evaluate for these 
examples aside of E/B, the strangeness per baryon s/B and entropy per 
baryon S/B as shown in the top section of the table 4. We do not enforce 
{s — s) = exactly, but since baryon asymmetry is small, strangeness is 
balanced to better than 2% in the parameter range considered. In the 
bottom portion of table 4, we present the compatible particle abundance 
ratios, computed according to the procedure developed in section 2. We 
have given, aside of the baryon and antibaryon relative yields, also the 
relative kaon yield, which is also well determined within our approach. 

The meaning of these results can be better appreciated when we assume 
in an example the central rapidity density of direct protons is dp/dy\ccnt. = 
25. In table 5, we present the resulting (anti) baryon abundances. We 
see that the net baryon density dh/dy ~ 16 it 3, there is baryon num- 
ber transparency. We see that (anti)hyperons are indeed more abundant 
than non-strange (anti) baryons. Taking into account the disintegration of 
strange baryons, we are finding a much greater number of observed protons 
dp/dy\'^^^ ~ 65 lb 5 in the central rapidity region. It is important when 
quoting results from table 5 to recall that: 

1) we have chosen arbitrarily the overall normalization in table 5 , only 
particle ratios were computed, and 

2) the rapidity baryon density relation to rapidity proton density is a 
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Table 4. For 7^ = 1.25, Ag = 1 and 7^, \q as shown: Top portion: strangeness per 
baryon s/B, energy per baryon E / B[GeM] and entropy per baryon S/B. Bottom 
portion: sample of hadron ratios expected at RHIC. 
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95 


110 


s/B 


18 
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0.87 
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0.16 
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0.16 
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0.15 
1.05 


A+A 
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Table 5. 
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21 


36 
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22 


11 


10 


0.7 


1.60 


1.03 


15 


25* 


21 


34 


30 24 


21 


10 


9.6 


0.6 



consequence of the assumed value of A^, which we chose to get E/B ~ 
100 GeV per participant. 

The most interesting result seen in table 5 , the hyperon-dominance of 
the baryon yields at RHIC, does not depend on detailed model hypothesis. 
We have explored another set of parameters in our first and preliminary 
report on this matter [49], finding this result. Another interesting property 
of the hadronizing hot RHIC matter as seen in table 4, is that strangeness 
yield per participant is expected to be 13-23 times greater than seen at 
present at SPS energies, where we have 0.75 strange quark pairs per baryon. 
As seen in table 5, the baryon rapidity density is in our examples similar to 
the proton rapidity density. 
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6. Conclusions 

We believe tliat the Fermi model interpretation of SPS strangeness re- 
sults decisively shows some interesting new physics. We see considerable 
convergence of the results around properties of suddenly hadronizing QGP. 
The key results we obtained are: 

1) same hadronization temperature r=142-144MeV for very different 
collision systems with different hadron spectra; 

2) QGP expected A^ = 1 for S and Pb collisions, and A^^ ~ 1.1 ; 

3) 7s > 1) indicating that high strangeness yield was reached before 
freeze-out; 

4) 7q > 1 as would be expected from high entropy phase and the asso- 
ciated value S/B ~ 40 ; 

5) yield of strangeness per baryon s/B ~ 0.75 just as predicted by gluon 
fusion in thermal QGP; 

6) transverse expansion velocity v^ = l/v3, the sound velocity of quark 
matter for Pb-Pb. 

Among other interesting results which also verify the consistency of our 
approach, we recall: 

• the exact balancing of strangeness {s — s) = in the symmetric Pb-Pb 
case; 

• increase of the baryochemical potential fiy^ = 203 it 5 > /i^ = 178 it 
5MeV as the collision system grows; 

• energy per baryon near to the value expected if energy and baryon number 
deposition in the fireball are similar. 

The universality of the physical properties at chemical freeze-out for S- and 
Pb-induced reactions points to a common nature of the primordial source 
of hadronic particles. The difference in spectra between the two systems 
arises in our analysis due to the difference in the collective surface explosion 
velocity, wf = 0.5 < v^ = l/vS, which for larger system is higher, having 
more time to develop. 

In our opinion, these results show that hadronic particles seen at CERN- 
SPS are emerging from a deconfined QGP phase of hadronic matter and do 
not undergo a re-equilibration after they have been produced. This has 
encouraged us to consider within the same computational scheme the pro- 
duction of strange hadrons at RHIC conditions, and we have shown that one 
can expect strangeness chemical equilibration in nuclear collisions at RHIC 
if the deconfined QGP is formed, with a probable overpopulation effect asso- 
ciated with the early strangeness abundance freeze-out before hadronization. 
We have shown also that (anti)hyperons dominate (anti)baryon abundance, 
and that rapidity distributions of (anti)protons are primarily deriving from 
decays of (anti)hyperons. 
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